Interplay between uniaxial strain and magnetophonon resonance in graphene 
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We study the fine structure of the phonon G peak in the inelastic light scattering spectra of 
strained graphene in the presence of a quantizing magnetic field. We show that under the conditions 
of the magnetophonon resonance (MPR), the spectral line shape in Raman spectra of the F-point 
longitudinal optic-transverse optical phonon undergoes multiple splitting and acquires a strong 
dependence on the polarization of the incoming and detected light. Interference of strain and 
MPR leads to the modification of the anticrossing branches, introducing the variation of the light 
polarization preferences depending on the proximity to the resonance, and their better resolution 
between neighboring anticrossings. 
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One of the most pronounced features in Raman spec- 
tra of graphene is a peak corresponding to the exci- 
tation of F-point longitudinal optic-transverse optical 
(LO-TO) phonons, known as a "G peak''.'^' Electron- 
phonon interaction affects the phonon dispersion and 
decaTjP"^, as well as couples these phonon modes to the 
lighlP. The renormalization of the LO-TO phonon by 
coupling with electrons is most pronounced in the pres- 
ence of a quantizing magnetic field which leads to the 
discrete spectrum of Landau levels (LLs) in graphene^. 
When energies of some inter-Landau-level transitions res- 
onate with the optical phonon energy, which is a condi- 
tion known as magnetophonon resonance, the electron- 
phonon (e-ph) interaction leads to the anticrossingline 
shape of the phonon mode in the Raman spectrurrPHil]. 
Besides the use of Raman spectroscopy for the analy- 
sis of the e-ph coupling in graphene, Raman measure- 
ments of phonon excitation can be used to probe strain 
in graphen^l^m^. Splitting of the G peak by strain has 
already been observed experimentalljW^"^, in agreement 
with density functio nal t heory (DFT) modeling of the 
Griineisen paramete j^° l ^^ l 

The renormalization of the optical phonon spectrum by 
the e-ph interaction is caused by a mutual transformation 
of the phonon into a virtual electron-hole pair and back. 
In a zero magnetic field this leads to the corrections to 
the dispersion relation of optical phonons and widening 
of the G pealP. In the presence of a strong magnetic field 
the spectrum of electrons is discrete, with Landau levels 
at e„± — ± j^\/2n {v is a Dirac velocity in graphene, 

and \b — \Jhcl eB is a magnetic length) in the con- 
duction (+) and valence (— ) band, and a peculiar zero- 
energy LL (n = 0) at the band edge. The discrete elec- 
tron spectrum makes coupling between the lattice mode 
and electron-hole excitations resonant when the optical 
phonon energy ujq — 1586 cm^^ = 197 meV matches the 
energy il„ = ^j2^(^\Jn -f 1 + ^Jn) of an e-h excitation 
with an orbital momentum An = ±1 and belonging to 
the representation E2 of the Cgi, symmetry group of the 



graphene latticJ^'^. As a result, magnetophonon res- 
onance (MPR) involves a pair of interband transitions 
(n + 1)+ (with An = +1) and (n -|- 1)~ n+ 
(with An — —1) with energies il^, and each of these 
two electronic modes is coupled to the phonons with op- 
posite circular polarizations.'^^] Here we stress that the 
MPR in graphene has a built-in natural selectivity of 
coupling with a circularly polarized (angular momentum 
An = ±1) optical lattice vibration, which is enhanced 
when one of the two resonant electronic modes is switched 
of by Pauli blocking in heavily doped graphene. In con- 
trast to MPR, strain transforms the degenerate F-point 
LO-TO optical phonon modes into a pair of phonons with 
energies Wg+ '^G- 1 linearly polarized along the two 
perpendicular axes of the strain tensor and split propor- 
tionally to the amount of strain, e. 



In this Rapid Communication we study the interplay 
between the fine structures in the optical phonon spec- 
trum induced by strain and the MPR, and predict its 
manifestation in magneto-Raman spectra of graphene. 
Since both of these effects are sensitive to the polarization 
of the lattice mode probed by Raman spectroscopy, the 
resulting MPR in strained graphene acquires a peculiar 
polarization and carrier density dependent form, which 
is analyzed in details for two major MPR conditions: F- 
point optical phonon resonance with — > 1^/1^ — > 
(at B - 25T) and 1" ^ 2+/2- ^ 1+ (at B - 4.5 T) 
inter-Landau-level transitions. 



The F-point optical phonon in graphene is character- 
ized by the relative displacement u = -^iP-A — ub) of 
atoms in sublattices A and B. Together, the LO-TO lat- 
tice modes, the electronic excitations and their coupling 
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are described using the Hamiltonian 



^n" ^n°^ : 



H — Hp\i + He + Hc-ph; He— 

na=0,l±,... 

HpY, = y (u^ +a;gu2) + 2Mujow ^ u„e„i.u^, (1) 
^^e-ph = -g\/2Mujo ^ (s„»,„3 • u)c+„ 



, Cm « . 



,j°,mPl=0,l±, 



Here, ujq is a "bare" phonon frequency, formally defined 
before the e-ph interaction is taken into account, and M 
is a reduced mass of carbon atoms. The phonon mode 
splitting by strain is determined by a strain tensor. 



■xy 



Rg' (o ^^]Re, Tre = 0, (2) 



Exx — £ COS 26, £xy = s sin 26, 



and angle 6 determines the rotation Re of the princi- 
pal axes of the strain tensor with respect to the crys- 
talline axes (zigzag direction) of the graphene lattice!^ 
The size of this splitting is characterized by the parame- 
ter w = ^d{ujG- —ujG+)/de « 13cm~Y% as estimated in 
Ref . 1181 The electron-phonon interaction is parametrized 
by the coupling constant g 0.2 eV^^. The selection rules 
for electron-phonon interaction in -ffe-ph are incorporated 
into matrix elements S(„_|_]^-|q„3 — (Sni^ (n+i)")* ^ such that 
S(„_|.i)Q„f) = for n > 1 and Sicg — ae(j, where 

®0/0 ~ ^(^^ ^^1*) unit vectors which determine 
circularly polarized modei^. 

To describe the renormalization of the phonon spec- 
trum and its manifestation in Raman inelastic scatter- 
ing, we estimate the full retarded Green's function of a 
phonon, 

/•OO 

D^,{uj) ^ 2muj^ (K(i),u,(0)])e'("+*)*di, (3) 
Jo 

where square and angle brackets denote the commutation 
with subsequent averaging over the ground state. Using 
a random phase approximation (RPA^ the latter can be 
described as 



,(0) _ 



73 7;? ""^ 



we. 



(4) 



Here, 7g stands for phonon relaxation rate, To ~ 2.7 eV 
is the A-B hopping parameter for graphene^, related 
within a tight-binding approximation to the electron ve- 
locity as 7o ~ ^hv/ac-c (oc-C is the distance between 
two neighbor carbon atoms in a graphene lattice), and 

^{g^. (e)G^, (e + c.) + (e)G«, [e + co)] (5) 



is a retarded polarization operator taken for given elec- 
tron density rig and magnetic field B which determines 
Pb — 2/7rAg density of electron states in a single LL 
(spin and valley degeneracies are taken into account). 
Electronic Green's functions are 



2/(^)1 [G« (6) 



where /(e) is a Fermi function. In principle, the sum 
over the electronic states in Eq. ([s]) is accumulated over 
the entire Brillouin zone. For the sake of the following 
analysis, this contribution can be separated into 

n^^{uj) = H„,.(o)lB^o„^=o + ^n,„'^(tj), 

where the value ^vi^{0)\b=o n =o "^^ polarization op- 
erator for undoped graphene in zero magnetic field at 
zero frequency can be used to determine the renormal- 
ized phonon frequencj*^ 



The latter can be calculated^^ using the tight-binding 
model, while all the features of low-energy electron states 
near the Dirac point in a magnetic field, doping, and the 
frequency dependence are incorporated in 



AH,, = (eo),(eo):AHo + (eo),(eo):AHo, 
AHo(w)= [AHo(-^)]* = 

2 -^r^ /(e(n+l)±) - /(CnT) 



(6a) 



q2 °° 

27r ^ 

n=0 



E 



UJ ± Qr 



+ 



S-^ ./(£(n+i)±) - .f(£«±) 



(6b) 



where fin = Vij^iVn^ 



1 + y'n) accounts for interband 



transitions, Q',^ — ^/2j^{y/n + 1 — y^) accounts for in- 
traband transitions, and spin and valley degeneracy is 
taken into account. The term corresponding to the first 
transition 1^ -o- comes with an effective factor 2, which 
is due to the difference of LL wave functions for n = 
and n > 1. 

Using group symmetry analysis, one can show that the 
optical r phonon, which belongs to the E2 representation 
of the Cgy symmetry group, couples to a pair of photons 
through the vectors d, determined in terms of polariza- 
tion vectors 1 (1) of an absorbed (emitted) photon and 
transfor ming according to the same irreducible represen- 
tation E^, 



v^i^—x.y \ ■^ y y/ 



Then, the spectral density of Raman scattering (normal- 
ized by the total Raman efficiency of the G peak) is de- 
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termined as 



/ = — Ini 
27r 



1 



detA 



J2 (de5)(de„)*A 



(8) 



det A = AcjcjA(^(^ — we, 

^00(00) =^G + AAn0(o) -uj~ ijG, 

-4oo = ^oo = w{exx + iexy). 

Using (de^j/Q) = ±i-\/2(leo/^)(le^/Q)*, we find that for 
circularly polarized right (O) and left (O) incoming and 
outgoing photons 



1 j-^ ^OO/OO 

TT detA ' 



O-s-O(O-i-O) 



0. 



(9) 



Using definitions in Eq. we find that for linearly po- 
larized photons, 



-Im 



^{Aoo+Aoo)+wecos{2ip+2(p+2e) 
detA '■ 



(10) 



where angles ip and ip determine the direction of the po- 
larization vector of incoming (1) and outgoing (1) pho- 
tons with respect to the crystalline axis corresponding 
to the zigzag direction in graphene. Finally, for non- 
polarized incoming light the total spectral density for the 
polarization-integrated Raman spectrum is 



1 Ago + A 

— im 

2tt det A 



The resulting LO-TO phonon spectral line shape 
formed by the interplay between the MPR (which pro- 
motes anticrossing of circularly polarized lattice modes 
with the inter-LL excitations) and strain (which pro- 
motes splitting of two linearly polarized phonons into in- 
dependent Raman signals) is shown in Figs. [I] and [2] The 
value of strain e = 1.2% used in these plots is such that 
the resulting phonon splitting w\e\ and the MPR anti- 
crossing energy g have comparable sizes. To illustrate the 
MPR-strain interplay, we plot the calculated grayscale 
maps of the phonon line intensity in polarization-resolved 
Raman spectra for in and out photons with circular and 
linear polarizations, which interrogate differently polar- 
ized lattice modes. The interplay between strain and 
MPR is most pronounced when the inter-LL transition 
(?i + 1)^ — > 71+ is forbidden by Pauli blocking of the 
Landau level while the transition to the Landau 
level (n + 1)+ is possible, which takes place in medium- 
doped graphene (middle column in Fig. [T]) . For doping 
rig = 1.5 X 10^^ cm~^ and magnetic field S w 27T which 
corresponds to the main ljq w \/2hv/XB MPR condition, 
the zero-energy LL is full, and the Landau level I"'' is only 
partially filled. Thus, in the vicinity of the anticrossing, 
only one of two circularly polarized phonon modes un- 
dergoes resonant mixing with the LL transition — > 1'^, 
whereas the other remains unperturbed (middle column 
in the top two rows in Fig. [ij. Away from the reso- 
nance, the phonon line is split into two linearly polarized 
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FIG. 1. The intensity of the Raman spectrum in graphene 
strained along the zigzag axis as a function of magnetic field 
B. These spectra are plotted for a phonon relaxation rate 
7g = 6cm~^ and strain e = 1.2%. Three columns correspond 
to different electron densities, and rows to different polariza- 
tions. Symbols O (O) stand for right- (left-) handed circular 
polarization of the incident and scattered light; || and _L de- 
note linear polarization of light parallel and perpendicular to 
the principal elongation axis of the strain tensor. 



branches, which is reflected by the linear-polarization- 
selective line shape in the bottom rows of Fig. [T] Note 
that for the density rig = 1.5x 10"'^^ cm~^, all anticrossings 
at lower magnetic fields are blocked by electrons filling 
completely the relevant LLs, whereas for the higher den- 
sity (right-hand column in Fig. [T]) even the MPR with 
the — l"*" transition is blocked. One may also notice 
a kink in the position of the phonon line in the latter 
case, resulting from a partial depletion of the level 1+, 
which activates the coupling of the — ?► l"*" mode with 
the phonon, downshifting the energy of the latter. Fi- 
nally, the MPR anticrossing in low-density or undoped 
graphene, shown in the left-hand column, displays an in- 
dependent anticrossing behavior of two linearly polarized 
phonons separately. This happens because couplings of 
phonons with — >■ 1+ and 1^ transitions are almost 
equal, and this degeneracy allows two linearly polarized 
phonons to mix with equal strength with two linearly 
polarized interband transitions. 

Figure [2] demonstrates the higher {n" — n+ -I- 1 and 
n~ +\ — , n ^ 1) LL anticrossings at smaller magnetic 
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FIG. 2. Low magnetic field spectrum for strain 
along the zigzag axis, phonon relaxation rate 7g 
Here, anticrossing 2~ — >■ 1"*" is blocked, and 1~ — ^ 
in Raman spectra. 



= 1.2% 



= 6 cm 
2^ is active 



fields. The abundance of those can only be seen in the 



lower-density structures, where the Fermi level for elec- 
trons is low enough, Ep < uic/^. The calculated spectra 
for this case are shown in the left-hand column of Fig. [2] 
where the MPR displays the independent anticrossing 
behavior of two split linearly polarized phonons, which 
is characteristic for the low-density structures. The size 
of strain is also important in determining the MPR line 
shape. For weak strain, w\e\ ^ g, the lineshape and 
density dependence of MPR is essentially the same as 
in non-strained structures analyzed in Ref. [8] For the 
dominating strain, w\e\ 3> g, each of the two linearly 
polarized phonons undergoes anticrossing with a linear 
combination of {n + 1)^ and (n -|- 1)^ — ;> n+ tran- 

sitions, resulting in the MPR fine structure displayed in 
Fig. [2] 

Besides the direct influence of strain on MPR via 
changing the properties of LO-TO phonons, one may con- 
sider the effect of static lattice deformations on MPR, 
via the modification of electronic properties of graphene 
in the nonlinear response to strairP^SSl and generation 
of a valley-asymmetric pseudo-magnetic field by inhomo- 
geneous lattice deformations^^. Out of those two, the 
influence of strain on Dirac velocity, making it slightly 
anisotropic, in linear approximation can be merely re- 
duced to the weak renormalization of the magnetic 
field. The role of inhomogeneous strain on the inter-LL 
transitions may be much more drastic, leading to inho- 
mogeneous broadening and an additional structure in the 
MPR line shape, as discussed in detail in Refs. iQlandlTTl 
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The rcnormalization of the Fermi velocity due to the strain 
in linear approximation is described by He = vav{Sv,y + 
C,£vv)Pu- New LL states have energy = -y/l — C,'^s^t^± 
and wave functions '4>'^± {x,Py) = ■\/xtp„± {xx, xpy), where 

^ = \J i+Ce ' ®° e-ph coupling and LL degener- 

acy remain unchanged. Thus, this effect for MPR can be 
reduced to the eflfective magnetic field B' = {1 — (^e^)B. 



